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Thin, solid bodies with metric symmetries admit a restricted form of reparameterization 
invariance. Their dynamical equilibria include motions with both rigid and flowing aspects. 

On such configurations, a quantity is conserved along the intrinsic coordinate corresponding 
to the symmetry. As an example of its utility, this conserved quantity is combined with 
linear and angular momentum currents to construct solutions for the equilibria of a rotating, 
flowing string, for which it is akin to Bernoulli’s constant. 

The mechanics of thin bodies are amenable to a description in which energies are effectively 
expressed with an integral over a surface or curve. On sufficiently long time scales, and when lon¬ 
gitudinal waves are not important, the effective action is often taken to be metrically constrained, 
with a low-dimensional body stress tensor enforcing local distances between moving pieces of ma¬ 
terial. This approach, either at the level of the action or equations of motion, has been used to 
describe systems at many scales, from the overdamped motions of flagella HUS, filaments ra , or 
macromolecules [5], to the inertial dynamics of elastic rods EHna, soft robot links nn , yarn JT2], 
flags HS, or sonar and telegraph cables in the ocean H2HE|. 

Thin bodies such as rods and solid membranes moving in a higher-dimensional space share 
characteristics with both rigid bodies and fluids. The metric constraints of inextensible solids 
preclude the broad reparameterization invariance inherent to fluids and fluid membranes, leaving 
only highly restricted reparameterization invariance associated with any metric symmetries the 
body may possess. Effectively one has global, rather than local, symmetries on the body. However, 
these intrinsic symmetries are not trivially identifiable with rigid transformations of embedding 
coordinates, because a lower-dimensional body can adopt extrinsic curvature in the surrounding 
space. Thus, unlike space-filling rigid continua, the body may display additional material flows 
in equilibrium in addition to, and distinct from, rigid motions. Simple examples include a string 
being unspooled [12], a cable being laid down on the ocean floor [14] . or the flowing skirt of a 
rotating dancer [32]. Figure [T] is a schematic of such a combined rigid-flowing configuration. 

We will describe such a system using several types of coordinates. Aside from a time coordinate, 
there are three important sets of spatial coordinates to consider. First, there are time-independent 
coordinates attached to pieces of material. For an inextensible body, these can also be associated 
with geometric parameters, such as arc length. Second, there are coordinates for which material 
motions tangent to the body have been removed. These are effectively Eulerian from the intrinsic 
viewpoint of the body, but they are attached to a body moving through an embedding space, so have 
some Lagrangian character as well. This viewpoint is called the “orthogonal gauge” in [18]. Finally, 
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FIG. 1. Two snapshots of a one-dimensional body X rotating with angular velocity ui z and carrying a 
uniform tangential flow of material described by the velocity field T9 S X, T constant. Shown are the material 
description X(s, t), which follows the flow on the body, and the body-Eulerian description X(?y, t), rj = s+Tt , 
which rotates rigidly. The two descriptions are identical if T = 0. Let the straight grey arrow on the left 
snapshot denote either description at some point in time t\. The snapshot on the right occurs at a later time 
f 2 ; the material position vector has moved along the curve, while both curve and body-Eulerian position 
vector have rigidly rotated. At this later time, the s value corresponding to the same r) value is smaller by 
an amount T(t 2 — fi); it corresponds to a piece of material that came from a location farther back on the 
curve. The snapshots correspond to the same range of 77 values, but different ranges of s values. 


there are the embedding coordinates, which we will assume to be the simple Cartesian coordinates 
associated with an inertial frame in E 3 . We will refer to these three sets of spatial coordinates as 
material or body (Lagrangian), body-Eulerian, and embedding (Eulerian), respectively. 

The action of a continuum may possess symmetries in any of these coordinates, and correspond¬ 
ing conserved quantities may be extracted by application of Noether’s theorems. Translational and 
rotational symmetries with respect to embedding coordinates correspond to familiar conservation 
laws for linear and angular momentum. For solids, symmetries with respect to material coordi¬ 
nates imply symmetries of the body metric. For space-filling inextensible solids PH, these must 
necessarily coincide with rigid motions of the embedding coordinates. Fluids have a greater range 
of allowable reparameterizations or “relabelings”, namely any transformations that conserve the 
density of labels (22ll25| . Some workers have investigated symmetries of thin solid-like or fluid¬ 
like bodies in a higher-dimensional background, including relativistic extended objects [26H29] , 
energy-bearing curves [30], or fluid membranes [3ll . 

In this note, we introduce a conserved quantity for thin body dynamical equilibria, correspond¬ 
ing to a symmetry with respect to either body or body-Eulerian coordinates. In the former case, 
which corresponds to rigid motions, the quantity is identical to a flux derived in a non-variational 
manner for uniform rods by Maddocks and Dichmann |32]|. In either case, it is a useful tool in 
the search for solutions. We illustrate these ideas with the concrete example of a rotating, flowing 
string, known as a “yarn balloon” in textile processing applications m- In this example, the 
conserved quantity is constructed from portions of the kinetic energy and body stress. Related 
problems involving axially moving structures such as belts have stimulated discussion of Eulerian 
action principles and the boundary conditions they produce [33j. Other related work seeks to 
extend action principles and Noether’s theorems to general non-material volumes [34h37| . 

Without loss of generality, we will discuss a one-dimensional body such as that of Figure |T[ Let 
us begin by describing the body with a Cartesian embedding vector X(s,i), parameterized by a 
material coordinate s and time t. Consider the action of a classical string, which is an integral over 
a constrained kinetic energy density 

C = lnd t X ■ d t X. - \a (<9. S X • S s X - 1) . 


(1) 
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Here, n is a uniform mass density, and the stress or tension a is a multiplier enforcing the metric 
constraint <9 S X • <9 S X = 1. The string is inextensible, and the material coordinate is also the arc 
length. The rather trivial metric of the body, <9 S X • d s X, is symmetric in s. 

The variation of body coordinates will be performed in an analogous manner to the time vari¬ 
ation that leads to the temporal boundary term representing the total energy (Hamiltonian) of a 
classical action. Recall that for an action of the form 

ft2 ft2 rs 2 

A= / dtL= dt ds£(X(s,t)) , (2) 

Jtl Jtl Jsi 

a time shift t —>• t + 5t results in a new action 


A + 5A = f d(t + 5t) f ds £ (X(s,t + St)) 

J t 1 J Si 


rt2—St 


rs2 


dt 


' t\—8t 


ds C (X(s, t + St)) , 


'Si 


= A-H‘~‘’St + 


r* i 

Jtl Js 


S2 xr 

ds — • d t X(s,t)St + O (6t 2 ) 

si 


(3) 

(4) 

(5) 


After integration by parts, the integral term may contribute to the temporal boundary term. In 
our case, 




- I 2 ds \^id t X ■ d t X. 
J Si 


'Si 

rs2 


ds\nd t X-d t X 


5t+f dt I ds fidtX ■ dt ( dtX6t ) — ad s X. ■ d s (dtX5t) , 

t=ti J S 1 

t=t2 pt2 / rS2 \ 

5t+ dtdt i ds ^dtX ■ dfXSt j + other terms . 

t=t\ Jtl \J Si J 


( 6 ) 

(7) 


We are only interested in the temporal boundary term, which has contributions from the shift in 
endpoints as well as an integration by parts of the kinetic term inside the integral^ This boundary 
term can be set equal to a constant, 


E = t “dslfidtX- d t X. 
J Si 


( 8 ) 


This is the total energy of the body, purely kinetic in this case because the constraints store no 
energy. 

As E is an integral over the entire body, it is not very useful in solving for the local dynamics 
of an extended object like a string. Instead, we would like to find quantities that are uniform 
along the string, that is, conserved along a body coordinate rather than in time. Now consider an 
analogous variation of the material coordinate s —> s + ds. We find 


[ t2 dt\-jr\ s = s s l5s+ [ 

Jt-\ L Js 


'ti 

pt2 

= - dt 

Jtx 


S2 5C 

ds — ■ d s X(s,t)5s 

si 

r s 2 

S=s 2 


lnd t X-d t X\ a a Z^ i Ss+ j dsd s (ad s X- d s X5s) 


'Si 


+ other terms. 


(9) 

( 10 ) 


If the system is symmetric in a spatial coordinate, the corresponding spatial boundary term will 
be a function of time that is uniform along the string. In the current context, that of dynamical 


1 The variation also yields bulk and spatial boundary terms that, for the present purely mechanical system, provide 
redundant information with respect to that obtained from variation of the embedding vector X. 
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equilibria, these functions of time will be proper constants, and will be referred to as such in what 
follows. For rigid body motion of X(s,t), a uniform reparameterization 5 s is a symmetry of the 
action, and yields the constant 


C s = UdtX -d t X + cr. 


(ID 


This is identical to the constant C 3 in [38], and to the constant \C in m if a is replaced with the 
appropriate component of the surface stress there. It is a particular manifestation of a conservation 
law for uniform, hyperelastic rods [32] , 

We may also vary the position vector, as if to derive momentum balance equations. After 
X —» X + 4X, we have 




rs 2 


5 A = / dt ds 


1 si 

rt2 r s 2 

/ dt ds d s (ad s X ■ <5X) + other terms . 
Jt\ J Si 


AC 

dX 


<sx. 


( 12 ) 


(13) 


If 5X corresponds to a continuous symmetry of the system, we may obtain new constants. For a 
system rigidly rotating about an axis z, there are two such symmetries, namely translations z5z 
along z and infinitesimal rotations (z x X)50 about z. Thus, 


P z = ad s X • z , 

J z = crd s X ■ (z x X) . 


(14) 

(15) 


These constants were called c\ and C 2 /oj in [38], and J z here is analogous to \ j z in m- 

We may now easily obtain a quadrature for rigidly rotating string equilibria. Note that we can 
form an orthonormal triad of z, 6 = |g*x|| 1 an( l r = ~z X 6, and that with the origin of X on the 
z axis, X = rr + zz, and, thus, z x X = r6 and 9 S X • r = d s r. For a string rigidly rotating with 
angular frequency u about z, (9*X = uirO, so a may be written as a function of a single variable 
r = ||z x X11. Now use the three constants 0, and ( |15[ ) to write: 

= 1 - (d s X ■ z) 2 - (d s X ■ Of , 

Pi + {PJrf 


(■ d 8 X • r; 

(< d s rf = 1 - 


( C s - 1 injj 2 r 2 y 


(16) 

(17) 


This can also be reformulated in terms of an elliptic integral in r 2 [38], but we need not do that 
here. Solving for r, one then obtains a directly from 0 . and the other two coordinates by 
quadratures using (14) and (15) or (16). Nondimensionalizing ( |17[ ), using some chosen length scale 
along with [i and w, will lead to a family of equations with our three constants P z , J z , and C s 
as parameters^] These three constants were found a century ago by Gray [39], who used them to 
construct a quadrature akin to ( |17[ ), but for the radial coordinate r in terms of the coordinate z 
along the axis of rotation, rather than in terms of the arc length s. 

Let us now consider rigid motion of a shape , with flow of material along it. In the case of an 
inextensible string, this must be a uniform flow along the tangent vectors, as illustrated in Figure 
[l]for a rotating shape. This is a trivial example of a Killing flow, a velocity field on a manifold that 
preserves its metric. Such rigid-flowing equilibria occur during the controlled, rapid unspooling of 


2 In a prior analysis of the rigidly rotating string [3S], one of these constants ( C s ) was absorbed into the others, 
effectively eliminating a scaling degree of freedom. In the analysis of a flowing skirt El, the shapes were presumed 
to have an additional symmetry, allowing description in terms of a spherical curve rather than a space curve; this 
yields equations with one less parameter. 
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yarn [12], Gray discussed properties of such a configuration, but did not construct the solutions 


The symmetry of this system is not in the material coordinate s, but in a body-Eulerian co¬ 
ordinate rj = s + Tt, where T is a constant tangential velocity on the string]^] We must consider 
a body-Eulerian action integral describing a rigidly moving shape into, along, and out of which 
material is flowing. We write 

rt2 rt2 rV2 

A= dtL r) = dt dr] C v (X(r/, t)) , (18) 

J t\ J t\ Jrji 

2 C n = n {d t X + Td v X) ■ (dtX + Td ri X) - a (d v X ■ d v X - 1) . (19) 

The velocity appearing in the kinetic energy term is made up of a piece dtX(j], t ) describing the 
motion— in our case, rigid rotation— of a point on the shape, plus a tangential flow Td^X{r], t) on 
the shape. If we were integrating over a material volume, the limits of integration of the r] integral 
would be time-dependent, and we would simply be rewriting the previous action. This is a different 
action, one that has a symmetry in ?/. The situation is reminiscent of many engineering problems 
involving a finite window of an axially moving system [33]. or a garden hose or other conduit m- 
Under the variation r/ —> r/ + dr], the spatial boundary term is 


Cr, + n (dtX + TdrjX) • Td v X - ad ri X • d v X , 
-\]id t X-d t X + \]xT 2 -a, 

-\nd t X • d t X - \]iT 2 -(a- ]xT 2 ) . 


( 20 ) 

( 21 ) 

( 22 ) 


The regrouping of terms between (21) and (22) anticipates subsequent expressions. It will be 


helpful to rewrite the current in terms of the original material description of X. Hence, 

C v - \]iT 2 = \]id t X{r], t) • d t X(r,, t) + a - ]iT 2 , (23) 

= \]i (d t X(s, t ) - Td a X(s, t)) • (5tX(s, t) - Td s X(s, t)) + a — ]iT 2 . (24) 


The constants C v and C s coincide when T = 0. Note the resemblance of 0 and ( |23[ ) to Bernoulli’s 
equation, applied to ideal fluid flowing along a fixed or rotating streamline. In fluid or solid, Coriolis 
terms arising from simultaneous flow and rotation do not contribute to the conserved quantities^] 
Because ]iT 2 is a constant, expressions (23) and (24) are consistent with the (unnumbered) ex¬ 
pression for the related constant C derived in a different manner in m- Several results in fluid 
mechanics, including conservation of potential vorticity and the Kelvin circulation theorem, are 
derivable either through relabeling symmetry or alternate approaches [221425] . 

The other constants are also modified by the flow, and our new choice of action integral. In 
terms of the body-Eulerian X(rj,t), 


P z = [~]IT ( d t X + Td, n X) + ad v x\ ■ z , (25) 

J z = [~]iT (d t X + Td n X) + ad v X\ • (z x X) , (26) 

and in terms of the material X(s,t), 

P z = (- fjTdtX + ad s X ) • z , (27) 

J z = (- fiTdtX + ad s X. ) • (z x X) . (28) 


3 This T is like the t of [T7], and should not be confused with the T that briefly appears in [3H] and there represents 
the total velocity in the direction of the tangents, including any contribution from rigid motion. 

4 Healey ED pointed out the equivalence of one-dimensional solids and fluids in his discussion of axially moving 
loops of string. 








6 


Equation (28) above corresponds with equation (21) in [IT]. Both expressions (27) and (28) will 
simplify further after explicit insertion of the material velocity c^X. 

The coordinate system constructed earlier is still useful here. We now have a shape rigidly 
rotating with angular frequency oj about z, with an additional uniform flow T along it. Thus, 
dt.X(s, t ) = urO + Td s X, and 


P z = (o- /jT 2 ) d s x. • z , (29) 

Jz = (cr — /uT 2 ) <9 S X • r6 — fiTuir 2 , (30) 

a — f.iT 2 = C v — \^T 2 — ijfjLu> 2 r 2 . (31) 


This new expression for J z already tells us that for nonzero T, generic solutions are three- 
dimensional, rather than coplanar with the axis of rotation. Note also the expression a — /LT 2 , 
reflecting the Routhian symmetry of inertial frames on sang. Without the non-inertial frame 
created by rotation cj, the equilibrium shapes would be unaffected by the flow T along the shape, 
the centripetal forces arising from this motion balanced merely by a shift in the tension a. Let us 
rescale using an arbitrary length scale L and the natural parameters of the problem: s —> s/L, 
r —> r/L , P z —> P z / /j,lo 2 L 2 , J z -» J z /fuu 2 L 3 , C v -> C' 7? ///w 2 L 2 , T —» T/ujL. Then, proceeding in a 
similar manner as before, we obtain 


(d s rf = 1 - 


P z + (Jz/r + Trf 


(32) 


Again, the stress and other coordinates follow directly or by quadratures. Equation (32) represents 


the full four-parameter family of equilibria of the rotating, flowing string. This new result extends 
prior work by several authors |38l oT)i EH1 - I46| j^] The geometries and stress distributions of the 
resulting curves are an interesting topic in their own right, and will be the focus of a future paper. 


The solutions (32) can also be obtained by manipulating variations of multiple auxiliary fields 
PH US], or through extensive calculations involving integrations of projections of equations of 
motion mm- However, here we have introduced a minimum of variables into the actions, and 
have not solved or even invoked an equation of motion. Symmetries and kinematical information 
are enough to directly reduce the problem to quadratures, using only a few lines of calculation. 
Presumably, this approach has general utility, and may be exploited for many other problems 
involving dynamical equilibria of thin bodies. Valuable extensions might include consideration of 
dissipative dynamics governed by nonconservative action principles m, or nonuniform potentials 
corresponding to the presence of a body force mmm- 
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5 Note that Mack’s [44] solution to this problem cannot be the full solution, because it reduces when T = 0 to 
straight lines perpendicular to the rotation axis. There are many other two-dimensional and three-dimensional 
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in a favorable light— it does not, in fact, coincide with the general coplanar limit in [38], but only with a single 
special case. Gray’s solution m was at that time unknown to the author and, apparently, to the others as well. 
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